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Abstract 

The direct CP violation parameter Re(e'/e) has been measured from the decay 
rates of neutral kaons into two pions using the NA48 detector at the CERN SPS. 
With 3.3 million Kl — > ir^n^ events collected during the 1998 and 1999 running 
periods, a result of Re(e'/e) = (15.0 it 2.7)xl0~^ has been obtained. The result 
combined with the published 1997 sample is Re(e'/e) = (15.3 ± 2.6)xl0~^. 
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1 Introduction 

The violation of CP symmetry was first reported in 1964 by J.H. Cliristenson, J.W. 
Cronin, V. Fitch and R. Turlay, who detected a clean signal of Ki vr+vr" decays 
CP conservation implies that the Ks and Kl particles are pure CP eigenstates and that 
Kl decays only into CP=— 1 and Ks into CP=+1 final states. The observed signal of the 
forbidden — > 27r decays (CP=+1) indicates that CP is not a conserved symmetry. 

CP violation can occur via the mixing of CP eigenstates, called indirect CP vio- 
lation, represented by the parameter e. CP violation can also occur in the decay process 
itself, through the interference of final states with different isospins. This is represented 
by the parameter e' and is called direct CP violation. L. Wolfenstein in 1964 proposed 
a super-weak force responsible for AS=2 transitions, so that all observed CP violation 
phenomena come from mixing and 6'=0. In 1973, Kobayashi and Maskawa proposed a 
matrix representation of the coupling between fermion families 0. In the case of three 
fermion generations, both direct and indirect CP violation are naturally accommodated 
in their model, via an irreducible phase. 

The parameters e and e' are related to the amplitude ratios 

_ A{Kl vr+TT-) _ , 

and 

"^'^ - A{Ks - vrOvrO) ' ' ' 

which represent the strength of the CP violating amplitude with respect to the CP con- 
serving one, in each mode. Experimentally, it is convenient to measure the double ratio 
R which is related to the ratio e'/e: 

R = TTTT^ FT-^TT/TTTTT ; ^ ~ 1 - 6 X Re[e e) (1) 

T {Ks ^ TV^TT^) T{Ks ^ TT+71-) 

By the mid-1970s, experiments had demonstrated that CP violation in the neutral kaon 
system is dominated by mixing, with the limit Kele'/e) < 10~^|Q. On the other hand, 
theoretical work showed that direct CP violation in the Standard Model could be large 
enough to be measurable This stimulated experimental effort with sophisticated de- 
tectors to measure Re (e'/e). The first evidence for the existence of a direct component 
of CP violation was published in 1988 In 1993, two experiments published their fi- 
nal results without a conclusive answer on the existence of this component. NA31 
measured Re{e'/e) = (23.0 ± 6.5)xl0-^ indicating a 3.5a effect. The result of E731 §, 
R.e{e'/e) = (7.4±5.9)xl0~^, was compatible with no effect. Recently, a new more precise 
generation of experiments announced results from samples of their total statistics. The 
KTeV collaboration measured an effect of Re{£'/e) = (28.0 ± 4.1)xlO~^ and NA48 
published a first result of Re(£:'/e) = (18.5 ± 7.3) x 10^^ ||10|- These observations confirmed 
the existence of a direct CP violation component. Current theoretical predictions are in 



the range up to ~30xl0 [|TI] 



This paper reports a measurement of Re(£:'/£:) with increased precision, performed 
by the NA48 experiment, using data samples recorded in 1998 and 1999. The correspond- 
ing statistics is seven times larger than that used for the published 1997 result 
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2 The method 

Measuring Iie{e'/6) to a precision of ~10~^ requires several million Kl and Kg —>■ 
TTTT decays. A sufficiently large flux of kaons is produced by the high intensity proton beam 
from the SPS accelerator. Data are accumulated using a fast and efficient data acquisition 
system, including a trigger with high rejection power and a large capacity data storage 
system. 

The design of the experiment and the analysis method focus on making the in- 
evitable systematic biases in the event counting symmetric between at least two of the 
four components of the double ratio. In this way, most of the important systematic ef- 
fects cancel to ffist order, and only the differences between two components need to be 
considered in detail in the analysis. This allows the systematic uncertainties to be kept 
sufficiently low. 

In order to exploit the cancellations, all four modes are collected at the same 
time and from the same decay volume. To achieve this, simultaneous Ks and Kl beams 
are produced in two targets situated at different distances from the decay volume. The 
intensity ratio is such that the tttt decay rates from the two beams are comparable. The 
kaon production angles are tuned to minimise the difference in the Kl tttt and Ks tttt 
decay spectra over a large range of kaon energies. The beam axes are almost coUinear, 
both pointing to the centre of the detector, so that the decay products illuminate the 
detector in a similar way. The similarity is further enhanced by weighting each Kl decay 
by a function of its proper time, such that the Kl decay distribution becomes almost 
identical to that of Ks- The small remaining differences in beam divergences and beam 
geometries are corrected using Monte Carlo simulation. The small difference between Ks 
and Kl slow beam intensity variations are eliminated by weighting the Ks events by the 
Kl/Ks intensity ratio. 

As a consequence of simultaneous data collection all losses associated with the 
detector, trigger and reconstruction efficiencies, and with the beam activity, cancel to a 
large extent between Kl and Ks- The ratio of tc^tt^ and tt+tt^ decays in each beam is 
independent of the absolute flux. 

Ks decays are distinguished from Kl decays by means of tagging. Protons directed 
to the Ks target pass through a high-rate tagging station (see section Ks events 
can be identified by comparing the registered proton time to the event time. Since this 
method is used for both the and the tt^tt samples, the double ratio depends only 
on the difference in the Ks misidentification probabilities between the two decay modes, 
and not on their absolute values. 

Backgrounds affect differently each of the four modes in the double ratio. High 
resolution detectors are employed to achieve an efficient background rejection. Small re- 
maining impurities due to three body Kl decays are carefully subtracted. 

After applying selection criteria, the four K ^ nir decay modes are counted in 
a common 70 to 170 GeV kaon energy interval. In this interval, the Ks and Kl decay 
spectra are similar to within ±15% (Fig. ^4|). In order to reduce the influence of the 
residual spectrum differences on the result, the events are separated into twenty bins of 
kaon energy, each 5 GeV wide. The event counts in both the tt^^tt" and tt^tt'^ modes are 
corrected for the ~10% probability of misassigning a Kl decay to the Ks beam due to 
the high proton rate in the tagging station. After applying all corrections bin by bin, the 
result is obtained by averaging the twenty double ratios. Remaining uncertainties from the 
reconstruction, mainly due to uncertainties in the detector calibration and the geometry, 
are evaluated using several independent methods. 
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3 Beams 




Figure 1: A schematic view of the beam hne (not to scale). 



The NA48 experiment is installed at CERN and uses the 450 GeV proton beam 



delivered by the SPS. Kl and Ks beams |12| are produced in different targets located 



126 m and 6 m upstream of the decay fiducial region, respectively (Fig. |1|). 

The SPS has a cycle time of 14.4 s. It is filled with protons in two batches, each 
11.5 /is long. The beam is accelerated to a momentum of 450 GeV/c with a 200 MHz 
radio-frequency system. It is then allowed to debunch, and is extracted by means of a 
slow excitation with 69 /xs periodicity (3 SPS revolutions) over a spill length of 2.38 s. 
The small remnants of the typical frequencies (200 MHz, 100 MHz, 87 kHz and harmonics 
of 50 Hz) surviving the filtering and debunching, together with burst-to-burst fluctuations 
of the extraction, result in an effective spill length of ~1.7 s. Since the Ks and Ki beams 
are produced concurrently, the Ks/K^ ratio is maintained stable throughout the burst 
to within ±10%. This ensures that both beams are nearly equal in their sensitivity to 
intensity variations of the proton beam. 



3.1 The Kl beam 

The primary high-flux proton beam (~1.5xl0^^ protons per pulse) impinges on 
a target (a 400 mm long, 2 mm diameter rod of beryllium), with an incidence angle of 
2.4 mrad relative to the Kl beam axis. The charged component of the outgoing particles 
is swept away by bending magnets. The neutral beam passes through three stages of 
collimation. The first "defining" collimator, placed 41 m after the target, limits the opening 
angle of the beam. It is followed, 63 m further downstream, by a second "cleaning" 
collimator, which prevents particles scattered or produced on the aperture of the defining 
collimator from reaching the detectors. The fiducial region starts at the exit of the "final" 
collimator, 126 m downstream of the target. At this point, the neutral beam is dominated 
by long-lived kaons, neutrons and photons. Only a small fraction of the most energetic of 
the short-lived component {Ks and A) survives. 
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3.2 The Kg beam 

The non-interacting protons from the target are directed onto a mechanicaUy 



bent mono-crystal of silicon |T3[. A small fraction of protons satisfy the conditions for 
channelling and are deflected following the crystalline planes. This component passes 
through a small aperture collimator incorporated in a beam dump, which absorbs the 
main, undefiected beam of protons. The fine-tuning of the beam focusing and of the 
crystal position and angle with respect to the incoming proton direction allows to select 
the flux of protons transmitted to be ~3xl0^ per pulse. Use of the crystal enables a 
deflection of 9.6 mrad to be obtained in only 6 cm length, corresponding to a bending 
power of 14.4 Tm. It acts only on the protons and guarantees a sharp emittance for the 
selected beam, without cancelling the deflection of the upstream sweeping magnet on 
other charged particles. 

After the beam dump-collimator, the transmitted protons pass through the tagging 
station (see section [4.1| ) which precisely registers their time of passage. They are then 
deflected back onto the Kl beam axis, transported through a series of quadrupoles and 
flnally directed to the Ks target (same size as Kl) located 72 mm above the beam axis. 
A combination of collimator and sweeping magnet deflnes a neutral beam at 4.2 mrad to 
the incoming protons. The decay spectrum of kaons at the exit of the collimator is similar 
to that in the Kl beam, with an average energy of 110 GeV (Fig. ^41) . 

The flducial region begins 6 m downstream of the Ks target, such that decays are 
dominated by short lived particles. At this point, the Ks and Kl beams emerge from the 
aperture of the flnal collimators into the common decay region. The whole Ks target and 
collimator system is aligned along an axis pointing to the centre of the detector 120 m 
away, such that the two beams intersect at this point with an angle of 0.6 mrad. The 
whole decay region is contained in a 90 m long evacuated tank with a residual pressure of 
~5xl0~^ mbar, terminated by a thin composite polyamide (Kevlar) window of 3x10"'^ Xq 
thickness. 



3.3 Particle rates in the decay region 

It is important to keep the neutral kaon beams as free as possible from the con- 
tamination of particles, such as muons and scattered neutrons and photons, produced in 
the target or along the beam line. The 2.4 mrad production angle is chosen to reduce the 
neutron flux per useful Kl decay. Muons are reduced by the use of the bent crystal and 
by a further order of magnitude by the subsequent sweeping magnets. Beam halo particles 
are suppressed by means of multiple collimation and veto counters. Only 2x10^ muons 
cross the detector for ~10^ Kl decays per spill, while the high fluxes of neutrons (~10^) 
and photons (~10^) in the neutral beam remain inside the beam pipe such that they do 
not reach the main detector. Since the fiducial decay region is enclosed in an evacuated 
tank interactions with air are minimised. 



4 Detectors 

The detector is shown in Fig. ^ [|I4|. Seven annular counters (AKL) surround the 
fiducial region to record photons escaping the acceptance of the main detector. Charged 
particles from decays are measured by a magnetic spectrometer. This is followed by a 
scintillator hodoscope which contributes to the trigger decision and also gives the precise 
time of charged decays. A liquid Krypton calorimeter is used to trigger and reconstruct 
27r° decays. It is also used, together with a subsequent iron-scintillator calorimeter to 
measure the total visible energy for triggering purposes. Finally, at the end of the beam 
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Figure 2: The detector. 



line, a series of muon counters are used to identify Kl—>^ tt/xz/ (K^s) decays. 
Two beam counters are used to measure the intensity of the beams. One is located at the 
extreme end of the Kl beam line {Kl monitor) and the other [Ks monitor) is at the Ks 
target station. 

4.1 The tagging station 

The tagging station (or Tagger) is located on the Ks proton path after the bent 
crystal. It consists of two scintillator ladders, crossing the beam horizontally and verti- 
cally [|15|. All scintillators are 4 mm thick in the beam direction and 15 mm long. To 
equalise the proton counting rate in all channels each ladder comprises 12 scintillators 
of variable widths, from 0.2 mm at the centre to 3.0 mm at the edges. An overlap of 
50 /zm between two successive counters ensures that the beam is completely covered. The 
scintillators of vertical and horizontal ladders alternate along the beam direction (Fig. 
A proton crosses at least two of them, one horizontal and one vertical. Photomultiplier 



17 18 19 20 



KS protons 



Figure 3: The arrangement of the counters of the tagging station. Counters 1-12 are 
oriented horizontally and counters 13-24 vertically. 



pulses are digitised by an 8-bit 960 MHz flash ADC module [TO. A ~100 ns window is 
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read out around the trigger time. The reconstructed time per counter has a resolution 
of ~140 ps, and two close pulses can be resolved down to 4-5 ns. Proton times are re- 
constructed offline by combining the information from horizontal and vertical counters. 
The coincidence of this time with the event time assigns the decay to the Ks beam. A 
demonstration of the principle is given in Fig. 4 for tt^tt' decays, where the difference of 
the event time from the closest proton time as a function of the vertex position in the 
vertical plane clearly shows the Ks and Kl beam assignment. 
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Figure 4: Time difference between event and its closest proton for Ks^ Kl tt+tt events 
as a function of the vertical position of their decay vertex. 



4.2 The magnetic spectrometer 

The spectrometer is housed in a tank filled with helium gas at atmospheric pres- 
sure. A thin evacuated beam tube allows the neutral beam to continue in vacuum. Two 
drift chambers (DCHl and DCH2) are located before, and two (DCH3 and DCH4) af- 
ter, the central dipole magnet. These chambers and their interconnecting beam tube are 
aligned along the bisector between the converging Ks and Kl beam axes. 

The integral of the magnetic field is 0.883 Tm, corresponding to an induced trans- 
verse momentum kick of 265 MeV/c in the horizontal plane. All components of the field 
have been measured. During the run, the current in the magnet is recorded and any rel- 
ative variation larger than 5xl0~^ is corrected for. The momentum scale is set adjusting 
the reconstructed invariant mass in tt~^tt~ decays to the nominal K^ mass. 

The drift chambers have an octagonal shape and an area of 4.5 m^ |jT^. Each is 
made up of four sets of two staggered sense wires planes oriented along four directions, 
each one rotated by 45° with respect to the previous one, and all orthogonal to the beam. 
This permits charged tracks to be reconstructed without ambiguities, and minimises the 
effect of wire inefflciencies by providing redundant information. The wire material of the 
eight planes and the filling gas correspond to 0.4% Xq per chamber. Only four planes of 
DCH3 are instrumented. 

With such large chambers, it is important to know the linear dimensions accurately 
and to have good control of the uniformity. The geometric accuracy due to the cumulative 
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uncertainty on the wire positions is better than 0.1 mm/m. In addition, the average plane 
efficiency is measured to be greater than 99%, radially uniform to ± 0.2% (Fig. ^.a). 
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Figure 5: Plane efficiencies for DCH1,DCH2 and DCH4 as a function of the track radius 
(a). Reconstructed 7r"'"7r~mass (b). 



A short drift distance (5 mm, corresponding to a drift time of 100 ns) contributes 
to the high rate capability of this detector. Track positions are reconstructed with a 
precision of 100 fim per view. The momentum resolution, obtained from the analysis of 
special electron runs, is <j{p)/p = 0.48%©0.009 xp%, where p is in GeV/c. The resolution 
of the invariant tt+tt" mass is 2.5 MeV/c^ (Fig. ^.b). The track reconstruction provides 
an event time with a precision of 0.7 ns. This is used to cross-check the main charged 
event time determined by the scintillator hodoscope. 



4.3 The scintillator hodoscope 

A scintillator hodoscope is placed downstream of the helium tank. It is composed 
of two planes segmented in horizontal and vertical strips and arranged in four quadrants. 
Fast logic combines the quadrant signals, which are used in the first level of the trigger 
for charged events. Offiine, a combination of hits located at the extrapolated tt+tt^ track 
positions is used to reconstruct the charged event time with a precision of ~150 ps. 



4.4 The liquid Krypton calorimeter 

The liquid Krypton calorimeter (LKr) is a quasi-homogeneous detector with an 
active volume of ~10 m^ of liquid Krypton. Cu-Be-Co ribbons of 40 fim x 18 mm x 
125 cm define ~13000 cells, in a structure of longitudinal projective towers pointing to 
the centre of the decay region []TS|. The calorimeter is 27 Xq long and fully contains 
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electro-magnetic showers with energies up to 100 GeV. The cross section of a cell is about 
2 cm X 2 cm, and consists of a central anode in between two cathodes. The electrodes are 
guided longitudinally through precisely machined holes in five spacer plates located every 
21 cm. They follow a ±48 mrad zig-zag in order to maintain the mechanical stability 
and to decrease the sensitivity of the energy resolution to the impact position. Good 
energy response is further guaranteed by the initial current readout technique which also 
provides a high rate capability. The signals are shaped to ~75 ns FWHM and are digitised 
asynchronously by a 40 MHz flash ADC |]19[. The dynamic range of the digitisers is 



increased by gain-switching amplifiers which change the amplification factor depending on 
the pulse height. The calorimetric information readout is restricted by a zero-suppressing 
hard-wired programmable algorithm to channels belonging to an energy dependent halo 
around the most energetic impact cells. The calorimeter was operated in a stable way at 
a high voltage of 3 kV. Around 0.3% of the cells were defective and are excluded from the 
analysis. 



The performance of the calorimeter |^ is studied using the electrons from the 



abundant Kl ireu (Kes) sample recorded along with the Re{e'/e) data taking. The 
electrons are used to improve the cell-to-cell uniformity from 0.4% after electronic cali- 
bration to 0.15%. To avoid possible correlations between non-linearity and non-uniformity, 
only electrons in the energy range 25-40 GeV were used for this purpose. Cell-to-cell cal- 
ibration factors were checked using photons from 7r° and t] decays in special runs. 

Fig. ^ shows the resolution of the ratio E/p of energy (from the LKr) and mo- 
mentum (from the spectrometer), obtained using Kes events after the inter-calibration 
procedure. Unfolding the measured contribution of the spectrometer to the momentum 
resolution, the following energy resolution is obtained, with E in GeV: 

a{E) (3.2 ±0.2)% (9 ± 1)9 



E Je E 



B (0.42 ± 0.05)% 

The 3.2%/V-E sampling term is dominated by the fluctuations of the shower fraction 
outside the cluster radius used in the reconstruction. The 1/E term is given by the total 
noise in the cluster. The constant term has contributions from the cell-to-cell calibration, 
the residual gap width variations, the response variation with impact distance from the 
electrodes and the pulse reconstruction accuracy. 

The energy response is linear to about 0.1% in the range 5-100 GeV. This is 
shown in Fig. 0, where the average E/p is shown for electrons. To account for losses in the 
material in front of the calorimeter, 45 MeV is added to the electron energy as computed 



by a GEANT-based Monte Carlo ||2^. This simulation also predicts a 0.05% non-linearity 
because of the slight change in the electrode gap width, due to the projective geometry 
of the calorimeter. This effect is only significant for high energy showers which develop 
deeper inside the calorimeter. 

The position resolution of the calorimeter is measured using monochromatic elec- 
tron beams. Comparing the electron centre of gravity position from the nine cells around 
the impact point with the extrapolated track point, one finds a resolution better than 
1 mm in both directions, for energies above 25 GeV (the average of the photon spec- 
trum). 

The calorimeter gives an accurate time signal for neutral events, which is then 
used together with the tagging station time measurements to distinguish Ks from Kl- 
The photon time resolution is of the order of 500 ps in the energy range from 3 to 
100 GeV. The 2Tr^ event time is known with a precision of ~220 ps. Tails coming from 
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Figure 6: Energy resolution of the LKr calorimeter, obtained with electrons from 
Kl — > nei' decays. 

misreconstructed times have been studied using a single Ks beam; they are below the 
level of 10"^ 



4.5 The Ks anti-counter 

The Ks anti-counter (AKS) is located at the exit of the Ks collimator. It is 
composed of a photon converter followed by three scintillator counters, and its main 
purpose is to veto all upstream decays from the short-lived beam. The converter consists 
of a 2.96 mm thick iridium crystal of 1.8 Xq PT|, which converts photons into electron 
pairs keeping the diffraction probability small. The crystal axis is aligned with the beam 
axis to maximise the pair production. Given the measured conversion efficiency, only 
2.1 xlO""^ of the accepted Ks — > 2tc^ sample are due to decays occuring upstream of the 
AKS. 

Charged particles, whether from photon conversion in the crystal or from tt^tt^ 
decays, are detected by the scintillator counters. The time is reconstructed offiine with 
a resolution of about 160 ps. The inefficiency of the measurement is due to the intrinsic 
counter inefficiency and the dead-time of the subsequent readout system, which might 
affect the detection of Ks — >■ 7r~^n~. This results in a total inefficiency of 1.0x10"^. The 
final double ratio is therefore corrected by (1.1±0.4) x 10~^ for the difference in the AKS 
inefficiency between tt^tt^ and tt+tt" modes. 

The AKS counter is used offiine to define the beginning of the decay region for the 
Ks beam. An event is rejected if a hit is recorded in time by the second scintillator. The 
AKS counter plays therefore an essential role in the experiment: the iridium crystal and 
the second scintillator give the geometrical references needed to control the energy scale 
and to define the fiducial region for Ks decays into tt^tt" and tt+tt^ respectively. 



4.6 Hadron calorimeter and muon counters 

A calorimeter made of iron and scintillators, 6.7 nuclear interaction lengths thick, 
measures hadronic energy. It contributes, together with the LKr calorimeter, to the total 
energy trigger. Following the hadron calorimeter are the muon counters. Three planes of 
scintillators shielded by 80 cm thick iron walls provide timing information that is used 
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Figure 7: The linearity of the LKr calorimeter, obtained with electrons from Kl —>■ neu 
decays. 

offline to identify the background from Kl -n-fiu (K^a) decays. 



5 Triggers and data acquisition 

The rate of particles reaching the detector is around 500 kHz. The trigger is de- 
signed to reduce this rate, with minimal loss from dead time and inefficiencies, to several 
kHz. A part of the read-out rate is reserved for redundant low-bias triggers that collect 
data used for the direct determination of the trigger inefficiencies. 

The trigger decisions are collected in a trigger supervisor system |^ that records 
the time of the accepted events, relative to the 40 MHz clock signal used to synchronise 
the experiment over its entire 250 m length [^. The combined trigger decision is sent 
back to all the read-out elements in the form of a time stamp. This time stamp arrives 
no later than 200 /is after the event occurrence, while all data are buffered in the front- 
end electronics. The time stamp is transformed into an address location in the buffer. 
A sufficiently large time interval (100-250 ns depending on the subdetector) around this 
location is read out and written to tape. In case of an occasional pile-up of positive trigger 
decisions, consecutive time stamps are queued to allow all subsystems to send their data 
to the event builder before receiving another time stamp. If, nevertheless, one of the 
subsystems fails to send the complete data, the whole trigger system is blocked, ensuring 
equal losses for vr+vr" and tc^tc^ decays. Event building is performed in an on-site PC farm 
and transmitted via a Gigabit network to the CERN computer centre, where another PC 
farm controls the tape writing and carries out the data reconstruction and the monitoring 
of the experiment. 



5.1 Trigger for tt^tt" decays 

The trigger for vr^vr^ decays p4| operates on the analogue sums of signals from 
2x8 cells of the LKr calorimeter, in both horizontal and vertical orientations. The signals 
are digitised, filtered to reduce noise and summed into 64 columns and 64 rows, thus 
providing two projections of the energy deposited. The summed energy E is computed 
along with the first and second moments of the energy distribution in each projection. 
Ml -r. Ml y, M2^x and M2^y. The moments are converted into kinematic quantities using a 
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"look-up table" system. The radial position of the centre of gravity, C — \JmI^^ + M^y/E, 
is an estimate of the distance between the kaon impact point at the calorimeter plane 
(had it not decayed) and the beam axis. The distance of the decay vertex from the 
calorimeter, D = {M2^x + M2^y)/E — C^/mx, where tuk is the kaon mass, is computed 
to determine the proper time of the kaon decay. In addition, the number of energy peaks, 
in space and in time, in both vertical and horizontal projections, is computed in bins of 
3 ns. 

The trigger requires an electro-magnetic energy deposit greater than 50 GeV, along 
with C < 15 cm and a decay vertex less than 5 Kg lifetimes (r^) from the beginning of 
the decay volume. Another requirement, that there are less than 6 peaks within 9 ns in 
both projections, helps to reject background from Kl — > Stt^. This condition is released, 
however, if accidental activity is detected close in time. 

The electronics is implemented in a pipeline which makes the trigger free of dead 
time. The resulting rate of this trigger component is 2 kHz with a latency of 3 /is. 

The efficiency is measured with events triggered by a scintillating fibre hodoscope 
placed at the depth of about 9.5 Xq near the shower maximum in the LKr calorimeter. It is 
(99.920±0.009)%, with no significant difference between Ks and Kl decays. Therefore no 
correction to the Re(e'/e) measurement is applied. The inefficiency is dominated by losses 
from an unresolved pile- up with accidental hits and by occasional energy mismeasurcmcnt. 
A cross-check of the efficiency measurement was performed with a very loose trigger 
condition {E >15 GeV) applied to events collected in single i^^-beam runs. 

5.2 Trigger for tt^tt" decays 

The 7r+7r~ decays are triggered with a two-level trigger system. At the first level, 
the rate is reduced to 100 kHz by a coincidence of three fast signals: 

1. Opposite quadrant coincidence in the scintillator hodoscope (Qx), where the quad- 
rants are defined with some overlap in order to avoid geometric inefficiencies. The 
remaining inefficiency (0.05%) is due to electronics and scintillator geometry, and 
is equal for Ks and Kl decays. However, a Qx signal cannot be produced in two 
consecutive 25 ns clock periods, which leads to a dead time of 0.5%. Signals in a 
sufficiently large time window arc recorded to allow the dead time to be applied 
(offline) on an event-by-event basis also to the vr'^Tr" sample. 

2. Hit multiplicity in DCHl integrated over 200 ns, requiring at least 3 wires hit in at 
least 3 views {2track). This has an inefficiency smaller than 10~^. 

3. Total calorimetric energy (Etot), made by summing the electro-magnetic energy from 
the 7r°7r° trigger with the hadron calorimeter energy, is required to be more than 
35 GeV. Owing to the low resolution in the time and size of the hadronic energy 
measurement, some good events fail to pass the energy threshold or the time co- 
incidence with the Qx and 2track signals. The Eto^ efficiency is (99.542 ± 0.018)% 
for Kl events (proper-time weighted) and (99.535 ± 0.011)% for Ks events. The 
correction for the efficiency difference is applied in bins of energy and amounts to 
(0.9 ± 2.2) X 10^^ on the average double ratio. 

A signal composed of Qx x 2track x Efot + Qx/D is sent to the second level trigger 
with a latency of 5 /is. The Qx/ D component, where D denotes a down-scaling factor, 
is added in order to measure the efficiencies of the 2track and Etot components. An 
additional down-scaled 2track x Etot trigger signal by-passes the second level trigger to 
allow efficiency measurements of the Qx signal and of the second level trigger. 
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The second level of the tt+tt trigger E^ consists of hardware coordinate builders 



and a farm of asynchronous microprocessors that reconstruct tracks using data from 
DCHl, 2 and 4. Triggers are selected if the tracks converge to within 5 cm, their opening 
angle is smaller than 15 mrad, the reconstructed proper decay time is smaller than 4.5 
Ts, and the reconstructed vrvr mass is larger than 0.95 rriK- The latency is variable but 
does not exceed 100 /xs and the output rate is 2 kHz. The efficiency is (98.319±0.038)% 
for Kl (proper-time weighted) and (98.353±0.022)% for Kg decays. The inefficiencies 
are due mainly to DCH wire inefficiencies (1.2%) with a contribution from algorithm 
imprecision (0.3%) and misreconstructions from accidental hits (0.2%). The correction for 
the second level trigger efficiency difference is applied in bins of energy with an average 
of (-4.5 ± 4.7) X 10"^ on R. 

The time available to extract the data from the spectrometer read-out ring buffers 
to the second level trigger is limited. This leads to a 1.1% dead time. The same dead time 
condition is applied to 7r°7r° candidates event by event, to ensure that the principle of 
collecting vr+vr" and 7r°7r° decays concurrently is respected. 

In order to avoid recording events with high hit multiplicity, an overflow condition 
is generated in the drift chamber readout whenever more than seven hits in a plane 
are detected within 100 ns. In this case the front end readout buffers of this plane are 
cleared and the time of the overflow is recorded. Both the second level trigger and the 
reconstruction of events are affected by this condition. 

Overflows are mainly due to showers induced by interactions of electrons or photons 
in the material surrounding the beam pipe in the region of the spectrometer and by 5-rays 
coming from interactions of charged particles with the drift chamber gas. Occasionally 
they are also generated by noisy ampliflers operated with low thresholds. In this case at 
most two neighbouring planes are affected. 

In the offline reconstruction, a window of ±312 ns around the event time is required 
to be free of overflows, both for tt+tt^ and ti^tt^ decays. This time window is larger than 
the sum of maximum drift time and reset time. In the tt^tt^ sample 21.5% of events 
are removed by the overflow condition which reduces the sensitivity to Ks/K^ intensity 
variations by an order of magnitude. 



5.3 Other triggers 

Several other triggers were collected continuously during data taking for systematic 

studies: 

— A trigger for 3tt^ decays, given by the down-scaled tt^tt^ trigger without the peak 
condition, used for Ks tagging studies. 

— A trigger for tt'^tt^) (where tt^ stands for the Dalitz decay vr^, ee'y) decays, com- 
bining the information from the LKr calorimeter and the spectrometer, used to test 
Ks tagging and the energy scale. 

— Beam monitor triggers used to record the accidental activity, with rates proportional 
to Kl and Ks decay rates. Beam monitor signals are down-scaled and delayed by 
69 fis which corresponds to the periodicity of the slow proton extraction (3 SPS 
revolutions). 

— Calibration triggers used to monitor and calibrate the LKr calorimeter, the tagging 
station and the scintillator hodoscope. 
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6 Data samples 

The NA48 experiment collected data for the Iie{e'/e) measurement in SPS running 
periods during the summer of three consecutive years: 1997-1999. The data from the first 
running period yielded 0.49 million — ^ 27r° events. The result was published in |l^ . 

In the year 1998, the total number of Ki — > 271° events collected during 135 days of 
running was 1.1 million. In the year 1999, an upgrade of the trigger and the event builder 
PC farm, as well as an increase in the operational stability of the detectors and electronics, 
contributed to obtaining smaller dead time and higher data taking efficiency. In addition, 
the SPS spill length was increased by 10%. This allowed the experiment to collect 2.2 
million Kl —>■ events in 128 days. In all run periods the polarity of the magnetic field 
was regularly inverted to allow systematic checks on vr+vr" decay reconstruction. 

In the year 2000, a special run took place to cross-check the Ks tagging systematics. 
In this run, the NA48 detector operated with vacuum in place of the spectrometer, and 
with the i^^-protons swept away after their passage through the tagging station. This 
allowed a direct measurement of the accidental coincidence rate between the Kl 27r°, 
37r° decays and the protons passing through the Tagger. 

Along with simultaneous Kl and Ks beam runs, several auxiliary runs were ded- 
icated to various systematic checks. Data with muon or with Kl beams without the 
spectrometer magnetic field were taken for alignment purposes. Scans with a monochro- 
matic coUimated electron beam were used for spectrometer and calorimeter calibration 
and alignment. Runs with either Ks or Kl beam only were taken on a regular basis to 
verify the Ks tagging performance and to check the data quality in single beam condi- 
tions. Tests of the LKr calorimeter calibration were carried out using data from runs with 
two thin polyethylene targets exposed to a tt~ beam producing 7r° —>■ 77, rj ^ 'j'-f and 
T] — >• 37r° decays at precisely defined vertex positions. 



7 Decay identification 

The raw data amount to 170 TBytes. After decoding and hit reconstruction, the 
data were filtered and compacted in several steps to allow for the iterative improvement 
of calibrations, alignment and corrections. The last step, before the final event selection, 
was sufficiently fast that it could be repeated several times following the refinements of 
the corrections. 



7.1 Reconstruction and selection of 7r°7r° events 

The reconstruction of ir^n^ events is based entirely on data from the LKr calorime- 
ter. The time and height of the pulses are measured using a digital filter technique. The 
first calibration is performed using a calibration pulser system. The cell response is inter- 
calibrated comparing energy and momentum (E/p) of electrons from Kes data. Further 
checks and fine tuning of the inter-calibration by ~0.1% were performed with vr^s pro- 
duced in a 7r~ beam by adjusting the reconstructed vertex of the photon pairs to the 
target position. 

Small drifts of the pedestal due to temperature effects are monitored and corrected. 
A pile up of signals within 3 /is causes the pedestal to shift. These shifts are detected by 
comparing ADC samples, stored in a buffer before the trigger, with the average pedestal 
level. If there is a significant difference between the two, then the stored samples are used; 
otherwise, the average pedestal is taken. This procedure minimises the influence of noise 
on the pulse height measurement. 
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Photon showers are found by looking for maxima in the digitised pulses from 
individual cells in both space and time, and accumulating the energy within a radius of 
11 cm. The shower position is derived from the centre of gravity of 3 x 3 central cells. Both 
energy and position measurements are corrected for a dependence on the distance of the 
impact point to the electrodes using data from electron beam scans. The transverse scale of 
the calorimeter is checked using Kgs-electron tracks and residuals of ~200 fim are applied 
as a correction to cluster positions. In order to account for deviations from the projectivity 
of the calorimeter, the cluster positions are recomputed at the shower maximum depth. 
The expected shower depth is estimated from Kes data comparing reconstructed and 
extrapolated electron shower positions and extrapolated to photons using Monte Carlo 
simulations. 

Overlapping showers are separated using expected shower shapes from a GEANT 
Monte Carlo simulation. The quality of the shower shapes is tested with data from electron 
beam scans. Based on these data, an additional correction is applied at the reconstructed 
shower level. Energy losses at the borders of the calorimeter are also accounted for, using 
information from Monte Carlo simulations and electron beam scans. Zero-suppression bias 
in showers with energy smaller than 5 GeV is reduced by parametrising the ratio of the 
energies deposited in a well defined 7x7 cell shower box and the total reconstructed 
shower energy, as a function of the photon energy and the number of cells read out. A 
small decrease of measured energy due to space charge accumulation during the spill, in 
average around 1.5 x 10""^, is corrected using Kes data The average energy loss of 
photons in the material before the calorimeter was determined with GEANT Monte Carlo 
to be 15 MeV. Shower energies are increased by this amount. A correction for residual 
nonlinearity of the energy measurement is derived from the parametrised E/p distribution 
of Ke3-electrons (Fig. 

Correct shower reconstruction is ensured by accepting only showers with energies 
between 3 and 100 GeV, and positions more than 15 cm from the calorimeter centre, 
11 cm from its outside borders and 2 cm from any defective cell. 

Any group of four showers, each reconstructed within 5 ns of their average time, 
is examined for the vr^vr^ signature. A minimum distance of 10 cm between the showers is 
required to resolve correctly the overlapping energy deposits. The event time is computed 
from an energy-weighted average of time measurements in individual cells, using the two 
cells with the largest energy deposits in each shower. Cells with signals more than 5 ns 
from the average time are excluded from the time computation. 

The distance d^ertex of the decay vertex from the LKr calorimeter is computed 
from the energies Ei and positions {xi,yi) of the four showers, with the assumption that 
they come from the decay of a particle with the kaon mass rrix moving along the beam 
axis: 

Je E EiE^ ((x, - x,)^ + {y, - y,Y) 

^ _ V ^ (n\ 

'-''vertex V J 

A resolution of ~50 cm is obtained on d^ertex- The vertex position z^ertex is equal to 
ZLKr — dvertex- The invariant mass of two photons TYlfyY IS then computed as 



ElE2((Xl-X2)2+ (7/1-1/2)2) 

"^77 = ^ (3) 

(^vertex 

The two m^^ masses (mi and 1712) are anti-correlated (Fig. H) because d^ertex contains 
information from all four showers. In order to find the best shower pairing, a variable 
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is constructed. In this variable, mi and m2 values are combined in order to remove the 
anti-correlation. 
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where 0"+ and (j_ are the corresponding resolutions, parametrised as functions of the 
energy of the least-energetic photon. The average values of 0"+ and (T_ are 0.42 MeV/c^ 
and 0.83 MeV/c^, respectively. 
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Figure 8: Distribution of Ks 27t^ candidates in the space of two reconstructed m 
masses, mi and m2. The contours correspond to increments by one standard deviation. 
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As electrons and photons have different energy losses, the final adjustment of the 
energy scale is done with photons. The energy scale is directly coupled to the distance scale 
(Eq. 0). It is adjusted by comparing the average vertex position of Ks 27r° candidates 
at the AKS edge with that produced by the Monte Carlo (Fig. with an accuracy of 
3 cm. The adjustment of the energy scale is applied as a function of the run period. 

The sensitivity of the double ratio to the energy scale is minimised by the choice of 
both the kaon energy range and the fiducial decay volume. At the chosen energy bound- 
aries, 70 to 170 GeV, the shape of the decay spectrum favours a cancellation of the losses 
and gains in case of an energy scale shift. 

The total number of selected decays is nearly independent of the distance scale 
because the acceptance is only weakly dependent on the longitudinal vertex position. For 
Ks decays the beginning of the decay volume is defined by the AKS counter. The end of 
the decay volume is chosen to be 3.5rs downstream of the AKS. A 3 cm shift in this cut, 
corresponding to the accuracy of the energy scale adjustment, together with the effects of 
energy boundaries and acceptance, would lead to an error of 2x10""^ in the decay rates. 

All systematic uncertainties in the tt^tt^ reconstruction are summarised in Tab. |l|. 
The residual deviation from linearity in the photon energy measurement was parametrised 
as: 

AE = a + (3E^ + ^rE (5) 

where E is the energy of the photon shower, and r is its distance from the centre of the 
calorimeter. An upper limit for each of the coefficients a, (3, 7 is determined on the basis 
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Figure 9: Distribution of the reconstructed vertex position Zy^rtex of the Ks —>■ tt^tt^ 
candidates at the beginning of the decay volume. The origin of the 
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the nominal AKS position. 



of the observed nonlinearities in Kes decays (Fig. |^), in iiT — > 27r°, 37r° decays, and 7r° and 
1] ~^ 11 decays from tt" beams. In addition, the stability of the energy scale is verified as 



a function of kaon energy (Fig. |T0|). The sensitivity of R to non-linearities is determined 
by modifying the reconstructed energies in the Monte Carlo. The modification takes into 
account also a non-linearity at < Q GeV observed in vr" —>■ 77 decays, not covered in 
the above parametrisation. 

The values of the shower depth used for the position adjustment are checked by 
comparing electrons from Ke3 with Monte Carlo simulations. The agreement is better than 
2 cm, which corresponds to a systematic effect of 1.2x10"^. The precision in the sharing 
of energies in overlapping showers is accounted for with a ±2.0 xlO"'' uncertainty on R 
which is the size of the correction from the electron beam scan. The transverse scale of 
the LKr calorimeter, relative to the spectrometer, was verified with Kes-electron tracks. 
The uncertainty on the double ratio of ±2.5x10"'^ corresponds to the difference between 
these data and measurement during the calorimeter construction corrected for the termal 
contraction. 

Due to hadron photoproduction in the liquid krypton, significant asymmetric non- 
Gaussian tails arise in the measurement of energies in about 3 x 10~^ of photons. These 
tails are parametrised using electrons and vr" decays and inserted into the Monte Carlo. 
The uncertainty associated to the effect of non-Gaussian tails is ±1.2x 10~^ on the double 
ratio. 

The scale adjustment is verified by studying 7r° and r] decays from a known origin 
produced by the vr" beam. Both the 7r° — >■ 77 and the 77 — 77 samples give consistent 
results. For the latter, the r] mass is determined from rj —>■ 3tt^ decays, cross-checking the 
method with Kl Sn^ decay^ The linearity of the distance scale along the beam axis 
is checked by comparing the results with the two polyethylene targets placed 1462 cm 
apart. The distance measured from data agrees within 1 cm with the nominal value. 
This is consistent with the uncertainties on the radial non-linearity 7, shower maximum 



To be published 
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Figure 10: Reconstructed AKS position as a function of kaon energy for Ks 27?*^ events. 



position and energy sharing from Tab. |l]. Another check is performed with K vr^vr^ 
decays, with subsequent vr" ee'j decay, comparing the vertex reconstructed from the 
e+e~ tracks to the one derived from the electro-magnetic showers. 



7.2 Reconstruction and selection of tt+tt" events 

The vr+TT^ events are reconstructed from tracks using hits in the drift chambers 
of the spectrometer. First, clusters of hits in each pair of staggered planes are found. 
The pattern recognition is based on hit positions without using the drift times. Clusters 
in DCHl and DCH2, within a range of 10-110 cm from the centre of the chamber, are 
assembled into track segments separately in each coordinate view (x, y, u, v). The segments 
are associated with space points in DCH4 with the constraint that the track is compatible 
with a straight line in the y view (vertical direction). 

Each track candidate is converted into one or more tracks using the drift time in- 
formation. Drift times are corrected for wire lengths and offsets using constants extracted 
from alignment runs with /i and beams without the magnetic field. Among tracks 



Table 1: Summary of systematic uncertainties on R in the ir^ir^ reconstruction. 





Units of 10"^ 


Energy scale 


±2.0 


Nonlinearity: 






a| < 1 X 10-2 GeV 


±2.5 




/3\<2x 10-^ GeV-i 


±2.3 




7 < 1 X 10"^ cm"^ 


±1.5 


at < 6 GeV 


±1.5 


Shower maximum 


±1.2 


Energy sharing 


±2.0 


Transverse scale 


±2.5 


Non-Gaussian tails 


±1.2 


Total 


±5.8 
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sharing the same space point in DCH4, the one with the best is chosen. After adding 
DCH3 information, the track momenta are calculated using a measured magnetic field 
map and alignment constants. 

A vertex position is calculated for each pair of tracks with opposite charge after 
correcting for a small residual magnetic field due to the magnetisation of the vacuum tank 
(~ 2x10^'^ Tm). The longitudinal vertex position resolution is about 50 cm, whereas the 
transverse resolution is around 2 mm. Only tracks with momenta greater than 10 GeV 
and not closer than 12 cm to the centre of each DCH are accepted. The separation of the 
two tracks at their closest approach is required to be less than 3 cm. The track positions, 
extrapolated downstream, are required to be within the acceptance of the LKr calorimeter 
and of the muon veto system, in order to allow for proper electron and muon identification. 

The time of the tt+tt^ decay is determined from hits in the scintillator hodoscope 
associated with the tracks. The events with insufficient information to determine the decay 
time accurately are discarded. This inefficiency is 0.1% and is equal for Ks and Kl- 

The kaon energy is computed from the opening angle 9 of the two tracks before 
the magnet and from the ratio of their momenta pi and p2 assuming a. K ^ tt^-k^ decay: 



The energy range and the decay volume are defined to be the same as in the vr'^vr^ mode. 
The small distance between the effective positions of the AKS converter and the AKS 
counter, (2.2 ± 0.1) cm is taken into account. 

In the vr+vr^ mode, the distance scale depends only on the knowledge of the relative 
geometry of DCHl and 2. Likewise, the uncertainty of the energy scale depends on the 
geometry, because the influence of the magnetic fleld uncertainty cancels in Eq. ^. The 
distance scale in the tt+tt" event reconstruction can be checked, as in tt^'tt^ decays, by 
comparing the position of the edge of the Ks tt+tt^ decay vertex distribution with the 
nominal position of the AKS counter. A small mismatch of 2.5 cm is found in the data 
corresponding to a shift in the distance of the two chambers of 2 mm or to a 20 /xm 
relative transverse scale mismatch. Neither of the two can be excluded. A correction of 
(2 ± 2)xl0~^, corresponding to a 2.5 cm distance scale change, is applied to the double 
ratio. 

7.3 Background rejection and subtraction in the tt^tt^ sample 

The background to the vrOvrO si gnal comes uniquely from Kl 3vr° decays. It 
is largely suppressed by requiring no additional showers within ±3 ns around the event 
time. In order to avoid possible losses due to noise, only showers with reconstructed energy 
higher than 1.5 GeV are considered for this cut. This does not affect the suppression of 
decays since, for kinematic reasons, very few photons with energy below 1.5 GeV coming 
from such decays can reach the LKr calorimeter. The remaining background consists of 
events with escaping or overlapping photons, resulting in only four reconstructed showers. 

Further suppression is achieved by requiring that the (Eq. is less than 13.5, 
which corresponds to 3. 7 a of the m^y resolution. The loss due to this cut is dominated 
by photon conversion and is ~7%. The residual background is evaluated by comparing 




(6) 



where 




(7) 
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the properties of the experimental data in the control region defined by 36 < < 135 
(Fig. |lT]), for both Ks and Kl decays, with the corresponding Monte Carlo distributions. 

In the Kg beam, the 4-shower events contain almost exclusively 7r°7r° decays, and 
the control region is populated by events which have an energy response or a shower 
coordinate displacement outside their normal Gaussian distribution. There are essentially 
two contributing processes. The first is a decay with a photon conversion (or a 7r° — ee7 
decay) where just 4 showers are reconstructed in the calorimeter. Since no cuts are imposed 
against the presence of charged particle tracks, such events can end up as tt^tt" candidates. 
Their value can be very large and, correlated with it, there can be a displacement of the 
reconstructed kaon decay vertex relative to the original longitudinal position. The second 
contribution is the occurrence, in about 0.3% of showers, of hadronic photoproduction in 
their early stage of development in the liquid Krypton, with a consequent reduction of their 
reconstructed energy. Both types of processes, leading to non-Gaussian tails, are taken into 
account in an appropriate Monte Carlo simulation. The simulation traces the electron- 
positron pair through the spectrometer and, for the effects of hadronic photoproduction, 
makes use of the experimental data obtained in calibration runs with electrons. The Monte 
Carlo for Ks gives a good fit to the experimental distribution over its full range, as 
well as to the distribution of the longitudinal decay vertex position for events in the 
control region. On the other hand, the experimental x^ distribution for clearly exceeds 
the one predicted by the (exclusively 27r°) Monte Carlo at high value of x^. 
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Figure 11: The comparison of the x^ distributions for and Ks 27r° candidates, 
showing the excess due to the 37?" background in the K^ sample. 

The background from 37?*^ decays populates, to a good approximation, the mi-m2 
space evenly, owing to its combinatorial character. The x^ is defined such that each bin 
corresponds to an equal area in the mi-m2 plane, which allows an almost fiat extrapolation 
into the signal region. An extrapolation factor of 1.2 ± 0.2 was derived from Monte Carlo 
simulation. This simulation, when normalised to the kaon fiux, agrees with the background 
level extracted from the data. 

The background is subtracted from the K^ —>■ 27r° sample in bins of kaon energy. 
The resulting correction on the double ratio, taking into account all uncertainties, is 
(-5.9±2.0)xl0-^ 
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7.4 

7.4-1 



sample 



Background rejection and subtraction in the tt+tt" 

Ks sample 

In order to eliminate the background from A — ^ pvr" in the Kg 
a cut is applied on the track momentum asymmetry, 



vr ' vr 



P-\ 



P++P- 



< min(0.62, 1.08 - 0.0052^^^) 



sample, 



(8) 



where Ek is the kaon energy in GeV. This cut is applied to both the Ks and the Kl 
samples, and its choice is motivated, apart from A rejection, by its improvement of the 
detector illumination symmetry. The residual A contamination was verified to be negligi- 
ble by comparing the invariant mpjr mass distributions of Ks —>■ tt^tt" candidates, with 
I^TTTT — between 3 and 5 sigma of the rrinn resolution, for opposite signs of the mo- 
mentum asymmetry. Due to the large A/A production asymmetry (7/1), in case of A 
contamination, the two samples would show a significant difference in the population of 
this control region which is not observed. 



TT+TT sample, the two semi-leptonic Kl decay modes, Kgs and K 



7.4-2 Kl sample 
In the Kl - 

are the dominant background sources. The Kes decays are suppressed by requiring E/p to 
be less than 0.8. K^3 decays are rejected when signals in the muon veto system, associated 
with the tracks, are found within 4 ns. Both cuts are also applied to Ks candidates in order 
to symmetrise losses of true tt+tt^ decays with a large electro-magnetic energy deposit 
(5.0%) or with a tt /iz/ decay (1.3%). 




Figure 12: Principle of the transverse momentum calculation viewed in the plane defined 
by the reconstructed kaon momentum px and the target T; V is the decay vertex and I is 
the reconstructed kaon impact point at DCHl. For a vertex misplaced due to resolution, 
Pt' remains independent of the target position. 



Further elimination of the semi-leptonic channels is achieved by requiring the in- 
variant mass m^TT to be compatible with the kaon mass to within 3a of the resolution, 
where the resolution is a function of the kaon energy. The cut on 771,^^ also rejects all direct 
emission nir'j decays [2S|, present only in the Kl beam. Most of the inner Bremsstrahlung 
decays are accepted; they are, however, Ks and Kl symmetric. This symmetry was tested 
by Monte Carlo simulation to hold to better than 10~^. 

The missing momentum carried away by the undetected neutrino in the semi- 
leptonic decays is reflected in pt- The variable pt is the transverse component of the 
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Figure 13: Comparison of the px'"^ tail of the — > tt+tt candidates with the sum of all 
known components. 



reconstructed kaon momentum with respect to its flight direction, assuming that it 
comes from the target. The Ks or Kl target is chosen according to the vertical vertex 
position. To equalise the effect of the resolution of the vertex position reconstruction for 
Ks and K^, the px' component of the kaon momentum (Fig. O) is chosen as a variable 



for background discrimination. This component is orthogonal to the line TI joining the 
target with the reconstructed kaon impact point at the DCHl. The relation between p^ 
and Pt is 

TV 

Pt = yjPt (9) 

where TV is the distance from the target to the decay vertex. By requiring p^'^ to be 
smaller than 200 MeV^/c^, most of the semi-leptonic background remaining after the 
m,r7r cut is rejected, while 'ir~^n~ event losses are small (0.1%). The asymmetry between 
Ks and Kl due to non-Gaussian tails is smaller than 2xl0~^. This is included in the 
reconstruction uncertainty in Tab. |^. 

In order to subtract the residual Kgs and K^3 background, two control regions are 
defined in the m,r7r-PT'^ plane. The first region, 9.5 < (m^^r — ^k) < 19.0 MeV/c^ and 
300 < pt'"^ < 2000 MeV^/c^, is dominated by Ke3 events, while the second, -17.0 < 
(m^TT — fTT-K) < —12.0 MeV/c^ and 300 < px''^ < 500 MeV^/c^, contains roughly equal 
numbers of Ke3 and K^3 events. Both regions are chosen such that they contain neither rni'-f 
events nor collimator scattered events (see section |7.5| ), and have sufficiently symmetric 
resolution tails in the Ks and Kl beams. 

To model the background distributions in the control and signal regions, a Kes 
sample is selected with E/p > 0.95, and a K^3 sample is obtained by reversing the muon 
veto requirement. The latter contains also tt ^ fiu decays which are taken into account by 
comparison with a similarly selected Ks sample. The number of true K^ vr+vr" decays 
in these regions is estimated from the Ks — > tt+tt" sample. The Ki tt+tt^ candidates 
are compared to the model samples and the scaling factors that best match the two 
background model samples are found. Their extrapolation into the signal region gives a 
background estimate of 10.1 xlO^'* for the Ke3 component, and 6.2x10^^ for K^3. The 
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independence of the double ratio on the choice of control regions has been tested, and all 
results are compatible within ±2xl0~^. 

As a further check, the p^'^ distribution of —>■ tt+tt" candidates is compared 
over a large pt'"^ interval with the sum of all contributing components (Fig. |13D, taking into 
account kaon decays from collimator scattering. Deviations of around 10% were shown to 
come from events with high hit multiplicity. Their influence on the background estimate 
is reflected in an increased uncertainty on the amount of subtracted background. 

The background subtraction is applied in bins of kaon energy and the overall 
correction on the double ratio is (16.9 ± 3.0)xl0~'^. 



7.5 Collimator scattering correction 

Both beams are surrounded by halos of particles from scattering in the collima- 
tors. Since the collimators are close to the decay region, the scattered particles manifest 
themselves through Ks —>■ tttt decays. 



7.5.1 Ks beam 

A beam halo in the Ks beam is formed by scattering in the collimator or in the 
AKS anti-counter. It is cut symmetrically in the tt^tt^ and tt+tt^ decay modes by requiring 
the radius of the centre of gravity (the distance of the virtual kaon impact from the beam 



axis), Cg = \jCg^ + Cgy, to be less than 10 cm. Cg is defined at the plane of the LKr 
calorimeter, where the two beams cross, and the cut is applied in all four modes. In the 
7r°7r° mode, Cg is computed from the positions (xj, yi) and energies Ei of the four showers. 

4 4 

E XiEi J2 ViEi 
r^oo _ 1 ^00 _ 1 /1n^ 

In the TT+TT" mode, the positions {xi,yi) are obtained by extrapolating the track tra- 
jectories before the spectrometer magnet to the plane of the LKr calorimeter. The Cg 
resolutions are almost the same for tt'^tt^ and n~^7r~ decays (Fig. 14 a). 



The value of the cut is chosen such that only a small part of the halo (~1% of 
all Ks events) is rejected. The acceptance differences between tt+tt^ and vr'^vr*' decays are 
taken into account in the Monte Carlo, which contains a parametrised shape of the beam 
halo. The accuracy of the cancellation is verified by varying the centre of gravity cut, and 
the uncertainty is taken into account in the acceptance correction. 



7.5.2 Kl beam 

The Kl beam is defined and cleaned by three collimators. Double-scattered events 
dominate the beam tails. Unlike in the Ks beam, in the Kl beam, the pr'^ cut rejects 
more scattered events than the centre of gravity cut. The pt'^ cut is applied only to vr+vr^ 
decays; in 7r°7r° decays, the scattered events are removed only by the centre of gravity cut 

(Fig. igb). 

The correction for this asymmetry is computed from reconstructed K^ — vr+vr" 
candidates with an inverted pt'"^ cut. The scattered events are identified by the rriT^T^ 
invariant mass, and the background continuum is subtracted (Fig ^). By extrapolating 
the computed kaon trajectory back to the planes of the final or cleaning collimator, rings 
having the collimator radii can be seen. This demonstrates the origin of these events. 
Scattered events have a 30% higher probability of containing an accompanying shower 
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Figure 14: Radial distributions of the centre of gravity for Ks and Kl beams after all 
cuts. 



not associated to the track than the decays of kaons coming from the target. Since such 
events are rejected in the selection, this factor has to be taken into account in the 
correction. The correction is applied in bins of energy and it amounts to — (9.6±2.0) x 10~^. 
The error is dominated by the uncertainty on the accompanying shower factor. 
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Figure 15: Collimator scattered events at pt > 200 MeV /c . 

One cross-check of this result involves searching for scattered events in the vr^vr^ 
sample with centres of gravity between 5 and 10 cm. This region is also populated by 
Svr^ background, and by decays from kaons scattered in the defining collimator, which are 
confined to low transverse momenta. A combination of the three contributions is fitted 
to the centre of gravity distribution of 271^ candidates. The model distributions of 

scattered events are obtained from vr+vr" decays, and the shape of the Stt" background is 
obtained from events with high x^, cross-checked with 5-shower events. 

A second cross-check consists of looking at the proper time distribution of events 
with centre of gravity between 5 and 10 cm (Fig. |TB|). The scattered events follow the 
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Figure 16: Direct determination of the number of collimator scattered events in the Ki — 
TT^TT^ sample with 5 < < 10 cm from the proper time distribution. 



falling Ks decay curve, whereas the 37r° background increases at the end of the decay 
volume. In addition, the distribution of events scattered at the defining collimator, which 
are a mixture of Kl and Kg decays, was simulated using the collimator geometries. 

All three methods give compatible results within the uncertainty of the correction. 



8 Ks tagging 

To associate an event with its parent beam, information from the tagging station, 
which records the protons that produce Ks particles, is used. A decay is labelled Ks if a 
coincidence is found between its event time and a proton time measured by the Tagger. 
The good time resolutions inherent to the detectors guarantee accurate identification. 
Fig. ITTI shows the tagging distributions for Ks and Kl decays to 7r"'"7r~ which have been 
identified as such by their vertex position in the vertical plane. A similar procedure is 
not possible for neutral decays. Tagging is therefore the only tool available to distinguish 
Ks from Kl. The selection of Ks and Kl samples is done by means of tagging for both 
the vr+TT^ and ti^ti^ modes, so that the uncertainties are kept symmetric. Two kinds of 
mistagging can occur: 

— A Ks decay can be assigned to the Kl beam. This probability is due to coincidence 
inefficiencies and it is denoted by asL- These inefficiencies are small (~10~'^), but 
they might differ for vr^vr" and vr^vr^ decays because the respective event times are 
reconstructed from different detectors. 

— A Kl decay can be identified as Ks- This is due to accidental coincidences between 
the event and proton times, and it is denoted as aLs- This probability only depends 
on the proton rate in the Tagger, so, to first order, it affects both 7r"^7r~ and tt^tt" 
decays equally. It is ~10% for the chosen coincidence window. 

Both mistagging fractions vary according to the width of the coincidence window. The 
tagging window used in this analysis, ±2 ns, was chosen to minimise the total uncertainty 
on R that comes from mistagging. 
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Figure 17: Time coincidence for charged Kg and Kl events, identified by their recon- 
structed vertex. 



8.1 Ks tagging inefficiency: cxsl 

The tagging inefficiency can be measured accurately in the charged mode, because 
the identification of the decay origin is possible from the vertex reconstruction. In identi- 
fied Ks decays, a fraction of (1.63 ± 0.03) x 10"^ lies outside the defined ±2 ns coincidence 
window, as shown in the distribution of the difference between the event time provided 
by the hodoscope and the closest proton time (Fig.|l^. Redundant time information pro- 
vided by the drift chambers demonstrates that ~80% of these tails are due to Tagger 
inefficiencies. The reconstruction inefficiency of the hodoscope time is responsible for the 
remaining ~3xl0~^. Since mistagging is predominantly due to the Tagger itself, it affects 
equally and TT+TT modes, so that its effect cancels in the double ratio. 

In order to measure the part of the inefficiency associated with the neutral time 
reconstruction, a large sample of Ks and Kl decays into 271^ and 37r'^ is used, where one of 
the photons converts into an electron-positron pair. The usual selection and quality criteria 
are applied, but without requiring the reconstruction of the complete decay. Accepted 
events are allowed to contain between four and seven clusters. The two tracks must cross 
to give a vertex and must fulfill the electron identification hypothesis, i.e. E/p > 0.95. 
The charged time is computed using the standard procedure applied to the two electron 
tracks. The neutral time is computed by the association of central and lateral cell times of 
all electro-magnetic showers satisfying some weak criteria which ensure that they belong 
to a single event. The event-by-event compatibility of the charged and neutral time is 
checked within the usual ±2 ns window (Fig. |T8|). 

The number of cases where the two time measurements differ by more than 2 ns 
corresponds to the differential inefficiency associated with hodoscope and calorimeter time 
reconstruction. Including the sensitivity of the neutral time computation to the number 
and to the energy of photons in the systematic uncertainty, one finds that the tails in 
the time estimator difference At are less than 2.5x10"^ at a 68% confidence level. The 
same data sample allows the precise determination of the absolute value of the neu- 
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Figure 18: Difference At between charged and neutral time estimators for neutral decays 
with a conversion. 

tral tagging inefficiency. Convoluting the neutral-charged time difference with the tt+tt"- 
— Tagger time difference, the tagging inefficiency for 7r°7r° is obtained. This inefficiency 
of (1.64±0.5) X 10^^, has been found to be in good agreement with a direct evaluation 
determined in single Ks beam, where all events are associated with a proton. A third 
method consists of checking the neutral tagging inefficiency in an identified Ks 2ti^ 
sample belonging to the main data set. For this, 7r°7r° data are used with a subsequent 
7r° — >■ ee7 Dalitz decay. The vertex position reconstructed from the electron and positron 
tracks allows Ks decays to be selected. All events outside the ± 2ns coincidence window 
correspond to Tagger or calorimeter reconstruction inefficiencies. The result is in good 
agreement with the two methods described previously, though it remains statistically lim- 
ited because of the small (1.2%) Dalitz branching ratio. Tab. ^ summarises the results on 
vr^vr^ tagging inefficiency. 

Table 2: Different evaluations of Ks tagging inefficiency for tt^tt^ events. 



Method 


a™ in 10-^ 


Conversions 


1.64±0.5(stat+syst) 


Single Ks beam 


1.90±0.42(stat) 


7r°7r° with Dalitz 


1.9^o:8(stat) 



The important consideration for the double ratio measurement is the mistagging 
difference between the charged and neutral modes. From the previous analyses, it can be 
concluded that there is no measurable effect within an uncertainty of ±0.5x10"^. 

8.2 Accidental tagging: olls 

Because of the high rate in the Tagger, K^ events can have an accidental coin- 
cidence with a proton. The probability for this to happen is proportional to the width 
of the tagging window. It can be measured directly in charged decays by looking at the 
fraction of those tt^tt" events identified from their vertex position as Kl decays, that 
also have a proton in the ±2 ns time window (Fig. p!9|a). On average, it was found that 
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a^g = (10.649 ± 0.008)%. Such a direct evaluation is not possible for neutral events. An 
indirect method is therefore applied to both the tt+tt^ and the vt^tt" modes, to evaluate 
the quantity AaLs = «is ~ (^ts ^ which leads to a correction to the double ratio. The Kl 
samples were identified as such by tagging, so that no proton is detected within ±2 ns of 
the event time. With this sample, the quantity Ptag is computed, which is the probability 
of having a proton within a ±2 ns sideband window, located before or after the event time 
(Fig. ITqId). These windows are centred at 5 ns intervals from the event time, and so follow 
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Figure 19: (a) The distribution of proton times with respect to event time for Kl tt+tt^ 
identified from vertex, (b) The variables ais, Ptag and W are schematically shown for 
Kl TTTT tagged events. 

the 200 MHz structure of the proton beam. Because /Stag is measured using the fraction 
of Kl events without an accidental coincidence in the central window, the proton rate is 
generally slightly lower than for the full sample. This is reflected in /Stag being smaller than 
Q^LS by a quantity W. Typically, W is ~10% of uls- Mainly due to the 100 MHz beam 
structure, /Stag varies with the window; however, the corresponding W value compensates 
for this, using the relation: 

- C^ts = i^^ag + W^') - {Ptag + W^-) 

The (3tag parameters are extracted from the same data that are used for the double ratio 
determination. In the 7r"^7r~ mode, can be computed by comparing and in 

all windows, using the Kl sample identified by vertex position. In a similar way, (3^^^ can 
be extracted. To measure W^^, the assumption is made that Kl 27r° and Kl 37r° 
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decays are tagged identically, giving 1^°° = W^^^. The abundant sample of Kl 3tt^ 
decays is then used to derive Finally, in order to increase the statistical significance 
of the measurements, several 4 ns wide windows are used in the sideband region. The 
results are shown in Tab. ^ 

Table 3: Measured values of parameters APtag, AW and Aa^s- 





Units of 10-^ 




3.0± l.O(stat) 


= - w+- 


1.3± l.l(stat) 


ActLs = afs - ^Ls 


4.3 ± 1.4(stat) ± l.O(syst) 



The quoted systematic error on Aols takes into account the small variation of the 
result depending on the set of out-of-time windows used in the computation. Several argu- 
ments support the assumption W^^ = W'^^^. Firstly, the trigger conditions for — > 27r 
and Kl —>■ 37r° are similar and highly efficient (~99.9%). Secondly, these two modes are 
reconstructed by the same detector with similar criteria. Finally, in 2000, a special run 
period was devoted to measure simultaneously both W^^ and W^^^. After passing through 
the Tagger, the Ks proton beam was deviated away from the trajectory leading to the 
Ks target, so that only particles from the Kl beam decayed in the fiducial region. All the 
tagging parameters analysed show good compatibility between the 2tt^ and 3tc^ modes. 
In particular it is found: 

VI/oo _ p^ooo ^ (^_g g ^ 2.9) X 10"^ 

which confirms the initial assumption. Even though there is an overall agreement of the 
tagging parameters between the 27c^ and 3tt^ modes, jS^^g, which is the largest part of af^g, 
is directly measured from the sample entering in the double ratio, namely 27r° events. The 
usage of the external Sir'^ sample is restricted to the determination of the parameter 
only. 

8. 2. 1 Origins of Aa ls ^ 

The measured difference Auls = (4.3 ± 1.8)xl0~^ (see Tab.^]) indicates that the 
7r°7r° and tt+tt" samples are recorded in conditions of slightly different intensities, although 
the same rejection of trigger dead time and drift chamber overflow is applied to both 
modes. The reason for this is the higher sensitivity to accidentals in the trigger conditions 
and reconstruction of tt+tt^ compared to tt^ti^. The flrst contribution to Aols, associated 
with 7i~^n~ events lost at the trigger, is (1.0 ± 0.3) x 10"'^. This has been verifled studying 
events collected with fully relaxed trigger. The second contribution, related to the event 
losses at the reconstruction, is studied superimposing good events from data or simulation 
to beam monitor triggers, recorded in proportion to the intensity (see section [10.1.1| ), and 
it is (2.5±0.4)xl0-^ 

Tab. 1^ shows the comparison between expected and measured AaLs values for the 
available data samples. It is noted that improvements in the tt+tt" trigger code and the 
lower instantaneous intensity allowed by a longer spill, result in a smaller loss for charged 
decays and consequently in a smaller AaLs in 1999 relative to 1998 data. The measured 
values of Aols and their variation between the two years agree well with the expectations 
within errors. 
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Table 4: Values of Aols in the two different data periods, and comparison with the 
expectations from known sources. 



Data Sample 


Expected Ack^s 
Units of 10-^ 


Measured Aa^s 
Units of 10-^ 


1998 


5.1±1.0 


8.3±2.9 


1999 


2.6±0.6 


2.6±2.1 


1998+1999 


3.5±0.5 


4.3±1.8 



8.3 Effect of the Ks tagging on R 

After all cuts, four event samples remain, namely tt+tt^ and tt^tt^ events tagged 
as Ks, and 7r°7r° and vr+vr" events tagged as Kl. The mistagged fraction of events is 
subtracted from Ks and added to Kl- The double ratio is then sensitive to AasL and 

ActLS- 

AR ^ -6{a'sl - a^,) AR ^ 2iafs - a^s) 

where the numerical factors depend on the absolute mistagging probability a^s and on 
the flux ratio Ks/K^. The corresponding corrections to the double ratio become: 

AR = (0.0 ± 3.0) X 10"^ for AasL 
AR = (8.3 ± 3.4) X 10"^ for Aols 

9 Acceptance correction 

Ks and Kl particles decay with very different lifetimes. Despite the use of a 
common decay region of 3.5ts, the decay distributions vary strongly along the beam 
direction. The resulting acceptance correction on R can reach about ±10%, depending on 
the kaon energy. To cancel the contribution from the different lifetimes to the acceptance, 
the Kl distributions are weighted by a factor W{t), where r is the proper time of the 
kaon: 

/(r from Ks target) 2it 

W[t) = 

/(r from Kl target) lix 

I{t) describes the complete kaon decay intensity into 2tt, containing Ks, Kl and inter- 
ference components: 

/(r) = e-^/^^ + Ir^pe"^/^^ 

+ 2|r/|Dpe-^^(i/^^+i/^s) cos(AMr - ^) 



For the kaon system decay parameters ts, tl, AM, rj and 0, the world averages 
given in the Particle Data Book 1^91 were used. The K^-K^ production asymmetry 



factor Dp is obtained by a fit to the data, in an energy range > 120 GeV where 
the interference term is significant. After weighting, the Kl and Ks decay distributions 
become nearly identical, as shown in Fig. The few Ks events reconstructed upstream 
of the AKS are due to resolution while in the Kl case the applied cut at ts=0 defines 
the beginning of fiducial region. A small difference in acceptances remains, related to the 
difference in Ks and Kl beam sizes and directions. This residual correction is computed 
using a large-statistics Monte Carlo simulation. 
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Figure 20: Reconstructed lifetime distribution for — > vr+vr events before and after 
weighting. 



By adopting the proper-time weighting method, the Kl statistics effectively de- 
crease and, within the fiducial region of S.Sr^, it results in a ~35% increase in the statis- 
tical error of R. 



9.1 Simulation tools 

The Monte Carlo precisely simulates the beams and detector apertures. The 
simulation of the beam includes the Ks component and the interference term, using 
the measured production asymmetry factor Dp. Particle interactions in the detector ma- 
terial are parametrised. The effect of kaon scattering and neutron interaction in the Ks 
collimator and the Ks anti-counter are included in the Monte Carlo by parametrising an 
energy dependent Ks beam halo, which reproduces well the experimental distributions of 
the centre of gravity for tt+tt^ and tt^tt^ modes. The measured wire inefficiencies are intro- 
duced before the event reconstruction. To describe the photon response in the calorimeter, 
a library containing 1.7x10^ showers is used, generated using GEANT in 36 energy bins 
from 2-101 GeV in order to guarantee good resolution and linearity and in a large 51x51 
cell area to allow realistic energy sharing for close particles. A parametrisation of the 
non-Gaussian tails is added to the photon energy, together with the measured electronic 
noise and known contributions from non-uniform response. Two other libraries are used: 
one for electrons to describe conversions or Dalitz decays, and another for charged pions. 

To obtain an accurate acceptance correction, 3x10^ kaon decays per mode were 
generated. After all cuts, statistics equivalent to five times the data were obtained. The 
generation is sub-divided into run periods that follow the spectrometer magnetic field 
direction, the beam position in the detectors, dead cells in the calorimeter and wire prob- 
lems, as they are observed in the data. Simulated events undergo standard reconstruction 
and analysis criteria. 



9.2 Correction on the double ratio and systematic checks 

The effect of the acceptance on the double ratio, as obtained from Monte Carlo, is 
shown in Fig. BT] in 5 GeV kaon energy bins. It is used to derive the acceptance correction. 
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For comparison, the acceptance effect without applying proper time weighting to the Kl 
events is also shown. 

The largest contribution to the correction comes from the difference between the Ks and 
Kl beams around the beam axis in the spectrometer for vr+vr" decays. The acceptance 
correction related to the n^ir^ mode is small. The illumination of the calorimeter for 
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Figure 21: Acceptance effect to the double ratio R as computed from Monte Carlo, in 
kaon energy bins. The case without weighting Kl data is also shown. 



ir^TT and tt^tt^ decays is shown in Fig. p2l in a comparison for Ks and weighted Kl 



events. The Ks/Kl ratio varies in the 7r~'"7r~ mode for track positions close to the beam 
tube (Fig. 0.b). This is reproduced by the simulation, so the effect is absorbed into the 
acceptance correction. In the tt^tt" mode, the photon distribution in the calorimeter is 
almost identical for Ks and weighted Kl events, as shown in Fig. |2^.d. 

The systematic uncertainties are evaluated by modifying certain parameters until 
the agreement between experimental and simulated distributions is affected. Varying the 
Ks beam halo, the K^-K^ production asymmetry, the beam shapes in the detectors and 
the wire inefficiencies, the correction is stable to within 3.3x 10~^. A systematic cross-check 
was done by sending the generated kinematic values of all 3x10® events per mode of the 
Monte Carlo through a GEANT based Monte Carlo in order to identify biases associated 
to the fast simulation. Within the statistical accuracy, no difference was observed in the 
neutral mode. In the charged mode, an event-by-event comparison showed a systematic 
difference of (-4.6 ± 2.3)xl0~^ 

While the effect is small, the average of the two evaluations, with a systematic 
error increased by 2.3x10"'*, is conservatively quoted as an acceptance correction. The 
final correction to R for the acceptance is: 

Ai?(acceptance) = (+26.7 ± 4.1(stat) ± 4.0(syst))xl0-^ 

10 Effects related to beam activity and detector response variations 

The identification of good K — >■ tttt decays might be affected by the presence of 
extra activity from various sources. Additional particles can either be produced simulta- 
neously with the kaon candidate in the target (called 'in-time'), or they can come from 
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Figure 22: The calorimeter illumination by tt+tt" tracks for Ks and Kl weighted events 
is shown in (a), while the Ks/K^ ratio for data and MC is shown in (b). Similarly, the 

calorimeter illumination for photons as a function of the radius for vr^'Tr'^ is shown in (c), 
with the Ks/ Kl ratio for data and MC shown in (d). Distributions are shown for events 
in the kaon energy range 95-100 GeV. 

uncorrelated kaon decays (called 'accidentals'). Accidental particles come mainly from 
the intense K^ beam, and the neutron and muon components accompanying it. In-time 
background is effectively suppressed for Ki decays given the beam collimation and the 
extensive shielding, while in the Ks beam a possible additional in-time component needs 
to be considered. Beam intensity variations do not affect the double ratio, thanks to the 
concurrent recording of all four modes. However, slow variations of the Ks/Kl intensity 
during the run that accidentally coincide with unaccounted variations in the charged-to- 
neutral event ratio (e.g. from a change in the detector response), would lead to a bias. 
The measured variance of the Ks/Kl intensity ratio is ~9%, while the 71^71^/71^11^ ra- 
tio, corrected for trigger inefficiency, varied by 1-2% during the 1998-1999 data period. 
The potentially resulting effect is taken into account in the analysis by weighting the Ks 
events by a factor extracted from Ks/K^ intensity variation monitoring. The effect of this 
weighting on R is small (< 3x10"'^), indicating that both the beams and the detectors 
were sufficiently stable during the two years of data taking. 

10.1 Accidental activity 

The pile-up of extra particles with a good event may result in an event loss, 
depending on the time and space separation in the detector. Owing to the simultaneous 
data taking, Ks and K^ events are equally sensitive to the accidental beam activity, to 
first order. A residual effect on R can be written as: 

AR = (Ai?)intensity + (^-R) geometry (H) 

The two terms have the following origins: 
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1. The intensity term, ( A i?) intensity, appears if there is a decorrelation between the 
two beams. Describing the Ks and Kl beam intensity time dependences as l5,L(t) = 
(Is,l) X (l+^s,L(t)), this would mean that ?7s(t) and rjiit) differ. If in addition the 
mean tt+tt" and tt^tt'' losses (A"*"" and respectively) are not equal an effect on 
the double ratio can be expected. Assuming that losses depend linearly on the beam 
intensity this effect is given by: 

(Ai?)intensity = (A^^ - A°°) X {^^^{^1 ' Vs)) (12) 

2. The illumination term, (Ai?)geometry, is non-zero if accidental losses depend on the 
detector illumination, since the latter is slightly different for the Ks and Ki beams. 
This can influence the double ratio if it is not equal for and vr^vr events. The 
corresponding effect can be expressed as: 

(Ai?)geometry = (5A+- - (5A°°) X (1 + {r]LVL)) (13) 

where 6X represents the difference in losses between the Ks and Kl beams. 
The measurement of event losses induced by extra activity is necessary for evaluating 
these terms. This is done through the overlay method. 

10.1.1 Overlay method 

Losses from accidental activity can be simulated by overlaying data with special 
beam monitor (BM) triggers. Beam monitors (see section |5.3| ) trigger in proportion to the 
Ks and Kl beam intensities, so that the BM triggers give a picture of the ambient activity 
and noise measured by the detectors at a given time. Detector information for the event 
and the BM trigger are superimposed, and the so-called overlaid event undergoes the usual 
reconstruction. For this purpose only BM triggers without DCH overflow condition are 
used. After applying the standard analysis procedure, the probability of losing or gaining 
an event is computed by comparing, event-by-event, the original and overlaid samples. 
This method allows the effect of the extra activity and noise that was present when 
the data was recorded to be reproduced and measured. However, the method depends 
on how well the BM trigger content describes the intensity variations seen by the data. 
A comparative study of the extra activity in BM triggers and kaon events showed an 
agreement that was better than 97%. 

To overcome the problem of doubling the detector noise in the overlaid event, 
leading to an overestimation of the losses, an alternative method is used in parallel: the 
overlaid Monte Carlo (OMC). This consists of using the BM triggers to overlay noiseless 
Monte Carlo events, instead of overlaying data. A small fraction of Monte Carlo events 
were overlaid consecutively with two different BM triggers (doubly overlaid Monte Carlo, 
or DOMC), in order to allow a direct comparison with the losses computed from overlaid 
data. 

With these methods, two sources of losses and gains for good events were identified: 
— The mortality due to a BM trigger. A BM trigger containing extra activity has an 
enhanced probability of killing an event in the overlay procedure. The extra activity 
depends on the beam intensity. Because of the different selection criteria, this would 
affect the tt+tt" and tt^tt^ modes asymmetrically. Examples of cuts sensitive to the 
intensity include the muon veto in tt+tt^ and the rejection of extra photon clusters 
in 7r°7r°. 
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— Noise in a detector. An event can be lost or gained after the overlay owing to migra- 
tion across cut boundaries. Certain selection criteria are sensitive to noise, like the 
distribution in tt'^tt^. Losses caused by noise do not depend on the beam intensity. 
The overall effect of the overlay is dominated by losses. The net effect of the overlay 
(losses minus gains) is shown in Tab. ^ OMC describes the reconstruction losses due 
to accidental activity and noise. Losses in the charged mode are dominated by the DCH 
multiplicity overflow condition, while in the neutral mode the cut accounts for the bulk 
of the losses. Results obtained by comparing DOMC and OMC are directly comparable 
to those from the overlaid data. The agreement is reasonable in the tt^tt^ mode. DOMC 
gives slightly fewer losses in the 7r"'"7r~ mode. This is associated with the smaller DCH 
multiplicity condition rate in DOMC. 

Table 5: Net effect of accidental activity, computed using data and Monte Carlo events 
overlaid with BM triggers. The LKr noise effect is subtracted. The DOMC column repre- 
sents differences between first and second overlay. 



Losses— Gains 


OMC 


DOMC 


Data overlay 


77 77 


0.8% 


1.5% 


1.4% 


77 + 77" 


1.8% 


2.5% 


2.8% 


~.-\-~.— —O—O 
77 77 — 77 77 


1.0% 


1.1% 


1.4% 



In order to explore the losses throughout the entire range of beam intensity, infor- 
mation from the beam monitor, placed at the end of the beam tube, is used. In Fig. 
the net accidental effect, computed from OMC, is shown for vr+vr" and events as 
a function of the recorded intensity at the time of the BM trigger used for the overlay. 
Losses increase faster in 7r"'"7r~ than in tt^'tt^ events. The overall behaviour is compatible 
with the assumption that the dependence of accidental losses with the rate is linear. 




0-2 0.4 0-6 0.8 1 1.2 1.4 1.6 



Intensity (arbritrary units) 

Figure 23: Losses as a function of the intensity seen by the BM trigger, computed by the 
overlay method. Lines are plotted to guide the eye. 

In what follows, the conservative value computed from data overlay is assigned 
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to the difference of the effect of accidental activity on tt+tt and vr tt events, namely 
X+--XOO= (1.4 ±0.7)%. 

10.1.2 Effect on the Double Ratio 

Having measured the accidental losses for tt+tt^ and tt^tt^ events, the two terms 
of Eq. [Tl| can be computed. The intensity variation term, (Ai?)intensity, depends upon 
how well the variations of the two beams {rjs, rji,) mimic one other, as well as upon the 
differential losses between the tt+tt^ and ti^tt'^ modes (A^^,A°°), according to Eq. 12. 



Several estimators have been used to study a possible difference in the activity of the two 
beams, based on accidental information in good events. Among these are the number of 
hits in the AKL counters, the number of accidental clusters in the LKr calorimeter, and 
the number of extra tracks, all of which have been measured in both Ks and events. 
The corresponding occurrence rates for these estimators are similar in the two beams. An 
independent method to measure the difference involves examining the correlation between 
the Ks and Ki beams. All results agree that, to within 1%, there is no intensity difference 
seen by Ks and Kl events. So from Eq. |12|we obtain (Ai?) intensity = ±3x10""^, where an 
additional, conservative factor of two is used in the uncertainty, to account for a possible 
non-linear dependence of losses with the intensity. 

The illumination term, (Ai?)geometry, is estimated by the overlay method with 
data or Monte Carlo events. Both computations agree that there is no significant ef- 
fect within available precision. The interpretation of the result is that all sources of 
losses are uniformly distributed in the acceptance. An upper bound on the illumination 
term can be quoted, according to the statistical significance of the measurement, namely 
(Ai?) geometry = ±3xl0~^. Finally, the uncertainty on the double ratio from accidental 
activity is: 

AR = ( A i?) intensity + ('^-R)geometry = ±4.2 XlO""^ 

10.2 In-time activity 

Particles produced in the Kl target, in time with the Kl are well suppressed by 
the 120 m of coUimation and magnets. In the Ks case, the 6 m collimator occasionally 
allows in-time particles to survive into the decay region, and to be detected together with 
a good kaon decay. 

To put bounds on such effects, the 7r°7r° sample accumulated in single Ks runs is 
used. The low intensity of the Ks beam ensures the absence of accidental particles. The 
7r°7r° mode was chosen for the clean identification of the number of in-time clusters; in 
the charged mode, the extended and irregular hadronic showers could bias the study. By 
removing the extra (in-time) cluster cut, the tt^tt^ signal is enriched with the following 
events: tt^tt^ events which have a photon conversion, giving 5-cluster final modes and 
7r°7r° events with additional in-time particles. Only in the case of photon conversion is 
the reconstructed deteriorated. By comparing the experimental distribution of 2tt^ 
events with five clusters to the one predicted by Monte Carlo solely for conversions, an 
excess of (0.7 ± 0.4)xl0~^ was found. This would correspond to the in-time activity, and 
has been cross-checked using various time windows for the extra-activity. Therefore, an 
upper bound on the in-time contamination in the Ks beam can be set to lxlO~^. This is 
included in the uncertainty on accidental activity in Tab. |^. 



35 



Table 6: Number of selected events after accounting for mistagging. 





Statistics in millions 




1998 


1999 


total 




1.047 


2.243 


3.290 


Ks ^ ttOttO 


1.638 


3.571 


5.209 


Kl — > TT + TT^ 


4.541 


9.912 


14.453 


Ks — > TT + TT^ 


6.910 


15.311 


22.221 


Statistical error 


18.0x10"^ 


12.2x10"^ 


10.1x10"^ 



11 Result 

The number of events passing all selection criteria, taking into account Ks mistag- 
ging, is shown in Tab. and their energy distribution in Fig. El. The final result is com- 



puted by dividing the data into 20 bins of kaon energy from 70 to 170 GeV, and calculating 
the double ratio for each bin using proper-time weighted Kl samples. All corrections are 
applied to each bin separately, and the results are averaged. In order to avoid a statisti- 
cal bias, instead of the simple weighted mean, the bins are averaged using a statistically 
unbiased logarithmic estimator. The average effects of the corrections on the raw double 
ratio, together with their uncertainties, are summarised in Tab. |^. 
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Figure 24: Energy distribution of selected events after accounting for mistagging and after 
proper-time weighting of Kl events. 

After applying all corrections to the double ratio, the result is 

R = (0.99098 ± O.OOlOl(stat) ± 0.00126(syst)) (14) 
The stability of the result as a function of various measurement parameters was studied 



extensively. The double ratio as a function of the kaon energy is shown in Fig. ^ The 
size of the assigned systematic uncertainties was tested by varying the most important 
selection cuts. By changing the background-rejection cuts, even large variations of the 
amount of subtracted background did not cause any shift in the result greater than the 
estimated uncertainty on the background. Similarly, no excessive deviations of the mea- 
surement were observed by varying cuts related to the acceptance, i^'^-proton tagging, 
uncorrelated beam activity, the decay volume and beam scattering (see Fig. EB). 
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Table 7: Corrections and systematic uncertainties on the double ratio. In order to obtain 
the effect on Iie{e'/e), the numbers must be divided by a factor of 6 (Eq. |l]). 
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Figure 25: Corrected double ratio as a function of the kaon energy. 



The stability of the result was also tested with respect to the time variation of the 
measurement conditions. In this context, the double ratio was measured in run periods 
defined by technical accelerator stops. Further checks were performed as functions of spill 
time, SPS revolution phase, 50 Hz mains phase, the spectrometer magnet polarity setting, 
and the time of day. 

Another test consisted of splitting the tt+tt^ data sample into converging and di- 
verging topologies, depending on the orientation of the track-charges in the magnetic field 
of the spectrometer, thereby selecting different acceptances and momentum resolutions. 
The steadiness of the drift-chamber read-out overflow rejection was tested by accepting 
a special class of overflow-events which affect only one pair of staggered drift chamber 
planes, and have only moderate effects on the trigger and reconstruction efficiencies. The 
variations observed in all tests are well within the assigned systematic uncertainties. 
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Figure 26: Stability of the double ratio with variations of the selection cuts. The grey 
band shows the uncertainty related to the cut concerned. -Rstandard represents the value 
from Eq. |14[ 

The data filtering, selection and calculation of the corrections were performed by 
several independent groups. One analysis adopting a different scheme of data compacting 
and filtering was performed in parallel to the one presented. In this analysis, several selec- 
tion criteria were chosen differently, and procedures (such as the background subtraction) 
were conceived independently. This analysis obtained a result fully confirming the above 
measurement. 

Converting the corrected double ratio using Eq. 0, the following measurement of 
the direct CP violation parameter He^e' /e) is obtained: 

Re{e'/e) = (15.0 ± 1.7(stat) ± 2.1(syst))xl0-^ (15) 

Adding the two uncertainties in quadrature gives: 

Re{e'/e) = (15.0 ± 2.7) x 10"^ (16) 

The combined NA48 Iie{e' /e) measurement, including the result from the 1997 data ||10|| , 
is: 

Re{e'/6) = (15.3 ± 2.6) x 10^^ (17) 
where correlated uncertainties are taken into account. 

12 Conclusions 

The NA48 measurement of the direct CP violation parameter R.e{e'/e), based on 
the analysis of data from the 1998 and 1999 data taking periods, is in good agreement 
with the result from the 1997 data. The combination of both results deviates from zero 
by 5.9 standard deviations. 
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